Capacitively coupled double quantum dot system in the Kondo regime 
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A detailed study of the low-temperature physics of an interacting double quantum dot system in 
a T-shape configuration is presented. Each quantum dot is modeled by a single Anderson impurity 
and we include an inter-dot electron-electron interaction to account for capacitive coupling that 
may arise due to the proximity of the quantum dots. By employing a numerical renormalization 
group approach to a multi-impurity Anderson model, we study the thermodynamical and transport 
properties of the system in and out of the Kondo regime. We find that the two-stage-Kondo effect 
reported in previous works is drastically affected by the inter-dot Coulomb repulsion. In particular, 
we find that the Kondo temperature for the second stage of the two-stage-Kondo effect increases 
exponentially with the inter-dot Coulomb repulsion, providing a possible path for its experimental 
observation. 

PACS numbers: 73.63.Kv, 72.10.Fk, 72.15.Qm, 73.23.Hk 






a 



o 



> 

O 



X 



I. INTRODUCTION 

Many-body electron-electron interaction is one of the 
most striking phenomena in low dimension condensed 
matter systems. In this context, quantum dots^ (QDs) 
have played a prominent role in the recent progress of the- 
oretical studies, ^'•^ as well as experimental realizations,^"* 
as they offer a unique opportunity for successful measure- 
ment of many-body-related physical phenomena arising 
at low temperature regimes. The relevance of electron- 
electron interactions in QD systems results from the 
strong confinement of the electrons due to the reduced 
sizes of typical structures. ^"^^ This interaction is respon- 
sible for several fascinating phenomena, e.g.. Coulomb 
blockade, ^"'^'"'^^ and Kondo effect, ^■-'^•^'^^ leading to char- 
acteristic behavior of the thermodynamical and trans- 
port properties, which depend drastically on the number 
of QDs, as well as on their topological configuration in 

the structure. In recent years, strong on-site interaction 
in double3'i3,i5-23 ^^^ ^j.jpie24-3i QD pQD j^^jj TQD) 

structures have received a great deal of attention when in 
the Kondo regime. However, on-site electron-electron in- 
teraction does not exhaust all the possibilities in multiple 
QD structures, as electrons can, due to their proximity, 
interact with each other, even when located in different 
QDs. It is important to note that recent advances in the 
lithography of lateral semiconductor QDs have allowed 
greater control over capacitively coupled DQD systems 
in parallel or in series. '^^"'^^ This makes it even more rele- 
vant to better theoretically understand the effects of ca- 
pacitive coupling over systems like the one studied in this 
paper. Only recently has this long-range interaction at- 
tracted more widespread attention of experimental and 
theory groups alike. ^^"^^ 

The T-shape configuration, where two QDs are mutu- 
ally coupled via a tunneling matrix element, while only 
one of them is coupled to metallic contacts, has attracted 
considerable interest as it allows the study of the two 



stage Kondo (TSK) effect, ^^'"'^"'^'^ by fine control of the 
inter-dot tunnel coupling. This effect results from the 
progressive screening of the localized spin of the elec- 
tron in each QD. Due to different effective couplings 
of the electron residing in each QD to the conduction 
band, these magnetic moments are screened at differ- 
ent temperature scales, which allows for the definition 
of two distinct Kondo temperatures. While the on-site 
Coulomb repulsion has received a great deal of atten- 
tion, the inter-dot Coulomb interaction in this particular 
system has not, to the best of our knowledge, been con- 
sidered yet. In this paper, we will explore how the TSK 
effect changes when a capacitive coupling is included be- 
tween the dots. We will use the Numerical Renormaliza- 
tion Group (NRG) method^^"^* to calculate the trans- 
port and thermodynamical properties of the DQD sys- 
tem. We will show that this new ingredient is responsi- 
ble for dramatic changes in the low-temperature physics 
of the system. In particular, as shown in Fig. 6, the 
Kondo temperature for the second stage of the TSK ef- 
fect increases exponentially with the inter-dot Coulomb 
repulsion. This may have important consequences to the 
experimental observability of this effect. The paper is di- 
vided as follows: In section II, we present the model and 
a brief description of the NRG method. In section III, we 
present numerical results for thermodynamical quantities 
(entropy and magnetic moment), and transport proper- 
ties (conductance). Finally, in section IV, we present our 
conclusions. 



II. THEORETICAL MODEL AND NUMERICAL 
METHODS 

We study a system composed by two QDs (from now on 
referred to as QDl and QD2) coupled by a tunneling ma- 
trix element as well as by capacitive inter-dot Coulomb 
repulsion. This system is described by the generalized 




FIG. 1: (color online) Schematic representation of the DQD 
system being analyzed in this work. Note the capacitive in- 
teraction, U' , between the dots. 



Anderson Hamiltonian, which can be written as 



H — Hi, 



Hcb + Hi 



hyh^ 



(1) 



where Himp is the Hamiltonian describing the QDs 
[which we define as "impurity region"], Hcb describes the 
conduction bands, and Hhyb describes the coupling of 
QDl and the conduction bands (see Fig. 1). More ex- 
plicitly. 



Hi, 



y^ £ic\„Ci^ + ^ UiUi^nii + U'nin2 



+t 



'EH 



4=1,2 



"{aC2a + CJ^aCla 



(2) 



where the operator c]^ (cia) creates (annihilates) an elec- 
tron in the i-th (i = 1,2) QD, with energy £,;, spin a^ 
Uia = c\aCia is the number operator, and ni = X^ct"'*'^- 
The second term in Himp corresponds to the on-site 
Coulomb repulsion, where, for simplicity we will take the 
intra-dot interactions Ui = U2 = U throughout this pa- 
per. The third term describes the inter-dot Coulomb 
repulsion U' due to the proximity of the dots, and the 
last term describes the coupling between the two dots, 
with tunneling matrix element t' . 



Hcb 



^^ikCfka^^ika, 



(3) 



NRG^^"^^ approach, which allows for a systematic assess- 
ment of the Kondo effect in impurity systems. Within the 
NRG, we logarithmically discretize the conduction band 
and map it into a tridiagonal form, which corresponds 
to a semi-infinite chain where the coupling between the 
sites has the form 



tjv 



(l + A-i)(l-A-^-i) 



2Vl-A-2^-Vl-A-2^-3 



^-N/2 



(5) 



where A is the discretization parameter (all results shown 
here are for A = 2.5). The "impurity" site {N — 1) pos- 
sesses sixteen degrees of freedom, corresponding to all 
the base-states necessary to fully describe the quantum 
state of the two QDs, while the other A^ > sites cor- 
respond to single non-interacting sites of Wilson's chain. 
Denoting the base-states of iJ_i = ffmip for the QDs as 

I 4>i)-i =1 mi,TO2)-i (6) 

where irik = 1,2,3,4 corresponds respectively to 0,t,l 
or tJ,. In this basis the Hamiltonian -ff-i has matrix 
elements 



]H. 



l\ij 



i7_ 



(7) 



By diagonalizing the matrix defined in Eq. 7 we obtain 
a set of sixteen eigenstates | ^i)-i with correspond- 



ing eigenenergies E, 



(-1) 



In terms of the base-states the 



eigenstates can be written as 



*»)-! 



T.^. 



(8) 



where Aij is the projection of the i-th eigenvector onto 
the j-th base-state. Once we obtain the eigenstates, we 
calculate all the necessary matrix elements for the next 
iteration, after which we add a new site A^ = 0. To 
describe the resulting system we enlarge the Hilbert space 
such that the new basis is constructed performing all 64 
possible combinations 



|0,)o :=|m)® I *j)_i. 



(9) 



where the operator cj^^ (c^fecr) creates (annihilates) an 
electron with momentum fc, energy eeki and spin a in the 
£-th lead {i = L,R). Finally, 
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hyb 
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k<7 
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(4) 



Notice that QD2 couples to the band indirectly through 
QDl. 

For simplicity, we assume the hybridization coupling 
V(,k = y to be real, independent of fc, and the same 
for both leads. The conduction band is characterized by 
a constant density of states given by Pc(w) = Q{D — 
\ijj\)/2D, where D is the half-bandwidth and <d{x) is the 
standard Hcaviside step function. To properly study the 
low-temperature physics of this setup we employ Wilson's 



where m = 1, • • • ,4 and j — 1,- ■ ■ , 16. This procedure is 
repeated until the system has reached its strong-coupling 
fix point. When the dimension of the Hilbert space be- 
comes larger than Ng, where typically A^^ = 2500, it is 
truncated by discarding the eigenstates corresponding to 
the largest eigenenergies. At each iteration A^ we keep 
the energy spectrum E^ , together with the matrix ele- 
ments necessary to calculate the relevant physical quan- 
tities. This procedure allows us to calculate thermody- 
namical properties, such as entropy S*, magnetic moment 
/i^ = ksTximp, spin-spin correlation (Si-S2), occupation 
numbers riic = {riia-), etc, as well as dynamical quanti- 
ties, like density of states (DOS) and conductance. For 
the entropy and magnetic moment it is usual to define the 
contribution from the impurity as Simp and fi 



quantities are generically written as Xi, 



imp- Ihese 
- X ~ Xq, 



where Xq is calculated in the absence of the impurity. 
Within the canonical ensemble, as a function of temper- 
ature T^v, we can write, 



X{Tn) = 



1 



■E 



Zn{Tn) 
where (3^ = (Tn)^^, and 
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is a characteristic temperature associated to the A^-th 
iteration, /3 is a real number of order 1, X is the operator 
associated to the quantity X and 



Zn{Tn) = y^e" 



■0nE 



(12) 



is the canonical partition function. 

The conductance is calculated by the generalized Lan- 
dauer formula 

/oo 
lm[%{^)][df{oj)/d^]dij, (13) 
-°° 

where %ioj) = 2ttV^ pciuj)G'(^iio) , Gq = {2e^/h) and 
Gii{u>) is the Fourier transform of the full interacting 
double-time Green's function 



GUt,t') ^ -zeit -t')([c,,it),clit')]) 



(14) 



which in the present case results to be spin independent. 
The Green's function is calculated at frequency ujn = 
£'(1 -h A"^)A"(^"^)/^ within NRG in a standard manner 
via Lehmann representation. 



G^ujn) = E 
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;* 
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(15) 



We also employ a logarithmic Gaussian broadening*^ of 
the discrete NRG spectrum in order to obtain a smooth 
curve for the QDs DOS at arbitrary frequency uj 



pliLo) = --ImGliLo), 



(16) 



necessary to calculate the conductance at finite temper- 
ature. 



III. NUMERICAL RESULTS 

In order to proceed with our numerical analysis, let us 
set D, typically the largest energy scale of the problem, 
as our energy unit {D = 1). We then choose for all 
calculations t/i = [/a == C/ = 5 x 10"* and F == 3.2 x 
10~^, so that U/ttTi « 5, where Fi — 2TrV'^po, and po — 



1/{2D). We will study in detail the two cases where 
[/' = and U' = [/, and also the range < U/U' < 1 
The bare levels £i and £2 will be controlled by the same 
gate voltage (Vg), such that ei = £2 = V^- With the 
parameters set above, and Vg = —U/2 (i.e., with the 
system at the particle- hole (p-h) symmetric point), we 
can estimate the Kondo temperature for the single QD 
as*7 Tk = ^/rv7cxp(-7^C//8Fl) w 2.89 x 10"'^ for QDl, 
when QD2 is completely disconnected (t' = and U' = 
0). 



A. f/' = case 

Although this case has been studied in great detail in 
Ref. 22, we will present below some results that will help 
us understand the more complicated situation at finite 
U'. In Figs. 2 and 3, we show results for the tempera- 
ture dependence of the entropy and the square of the to- 
tal magnetic moment of the DQD system, Simp/ks ln(2) 
and ^iirnp/{9^^BY, for U' = 0, £1 = £2 = —U/2 and var- 
ious values of i', where /cg is Boltzmann's constant and 
p^B is the Bohr magneton. For the special case where 
t' ~ [(black) Ocurve], the DQD corresponds to the 
case where just QDl is coupled to the conduction band 
and QD2 is completely decoupled from the rest of the 
system. In this situation, as the temperature decreases, 
we find the following regimes: (i) ior ksT ^ D ^ U the 
DQD is in its free orbital (FO) regime; in this regime, 
the temperature is high enough to allow for all the six- 
teen DQD states to be populated. This results in an 
entropy Simp = fcs ln(16) (see Fig. 2) and a total square 
magnetic moment p? = 2 x [gpsY 1^ (Fig- 3), (ii) for 
ksTii < ksT < £i,£2, thermally excited charge fluc- 
tuations are suppressed, and therefore the entropy de- 
creases to Simp = fc_B ln(4) , as only states with one elec- 
tron in each QD are favored. This implies that p^ in- 
creases, as double and unoccupied states in each QD are 
suppressed, and the DQD is in the so-called local mo- 
ment (LM) regime. As the temperature decreases fur- 
ther, and becomes lower than Tk, QD2 remains in its 
LM regime, while the other spin is progressively screened 
by the conduction electrons due to the formation of the 
Kondo state, with the electrons in the leads screening 
the spin of QDl. This regime is characterized by the 
plateaus S^mp = fcBln(2) and /xf,„p = {gpsY/^. These 
contributions arise just from the spin in QD2, as it is 
never Kondo screened when t' = 0. 

For finite t' , however, the behavior of the DQD for 
temperatures below Tk is quite different from the one 
just discussed above for f' = 0. For very small values 
of <', as shown in Ref. 22, the DQD system presents a 
TSK effect, where both QDl and QD2 spins are screened 
by the conduction band, but the screening of the spin in 
QD2 occurs at a much lower temperature than that at 
which the spin in QDl is screened. The second screening 
stage emerges at a characteristic temperature T'j^ which 
depends mainly upon t' and Tk (the Kondo temperature 
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FIG. 2: (color online) Temperature variation of the entropy 
for various values of t' (see legend). The other parameters are 
[/i = f/2 = 5 X 10"*, U' = 0, ei = £2 = -2.5 x 10"*, and 
V — 3.2 X 10""^. Note that the temperature in the horizontal 
axis is scaled by the Kondo temperature of QDl (see text), 
which, for the parameters chosen, is Tk = 2.89x 10~^. Details 
for the different regimes are provided in the text. 
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FIG. 3: (color online) Magnetic moment ksTx as function of 
temperature for various values of i'. All parameters are the 
same as in Fig. 2. 



for QDl) as 



T'k 



aTKC:xip{-bTK/J'), 



where a and b are real positive numbers, with values 
approximately 1, and J' = 4i'^/[/ is an effective anti- 
ferromagnetic coupling (between electrons in QDl and 
QD2) favoring a local singlet state that competes with 
the regular Kondo energy scale Tr-. As a result, a TSK 
effect is expected for J' < Tk- We can clearly see this 



behavior, for example, when t' /U = 0.0030 [(blue) O 
curve] in Figs. 2 and 3. Note that for this value oi t' /U 
the second drop in the entropy (signaling the screen- 
ing of the electron in QD2) happens at T^ « 10~^^, 
with J' « 1.35 x 10~^ (for this value of t') being much 
larger than Tj^ (remember that Tk ^ 2.89 10"'^. As 
t' increases, and eventually J' becomes larger than Tk, 
a strong singlet is formed locally, destroying the TSK 
picture. A detailed discussion of this crossover can be 
found in Ref. 22, here, instead, we focus on the effect of 
the inter-dot Coulomb repulsion, as discussed in the next 
section. 



B. Finite U' 



1. Entropy and magnetic moment 
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FIG. 4: (color online) Same as in Fig. 2, but now for U' — U 
and different values of t' (see legend). The Kondo temper- 
ature used to scale the horizontal axis is that obtained for 
U' = U and i' = 0. The inset contains a comparison of f' = 
results between U' — [(black) O curve] and U' = U [(red) O 
curve] , showing that the Kondo temperature of the first stage 
(for t' = 0) does not depend on U' . 

Now, we turn on U' and for simplicity we choose U' = 
U (later on, we will analyze results for < U' < U). The 
new features to note in the entropy (in relation to the 
results in Fig. 2, for U' = 0) are (i) the very clear plateau 
(in Fig. 2) at Simp/kB ln(2) = 2 splits into two narrower 
(17) plateaus in Fig. 4: one at Simp/kB lii(2) w 2.66, and 



the one at Simp/kB ln(2) = 2 becomes less well defined, 
more like a shoulder. The plateau at sa ln(6)/ln(2) — 
2.58 (starting at T « 10^ Tk) comes from the fact that 
now, as U' = U, all 6 states with 2 electrons have very 
similar energies. As the temperature goes further down, 
the 4 states that can participate in a Kondo state in 
QDl (one electron in each QD) will have lower energy 



and the shoulder around Simp/kB ln(2) = 2 will form (at 
T m Tfc), (a) to obtain the second stage of the TSK 
(signaled by a plateau at Simp/kB ln(2) = 1, for finite t') 
one needs to go to one order of magnitude smaller values 
of i', when compared to the results for [/' = (compare 
the (purple) A curve in Fig. 2, for t' = 0.004, with the 
(purple) A curve in Fig. 4, for t' = 0.005, where all traces 
of the second stage in the TSK effect have vanished), 
(in) it is also apparent that, contrary to the t/' = 
case, Tjj. depends much more strongly on t' . Indeed, 
it is clear that the temperature at which the entropy 
starts to decrease to zero becomes considerably higher as 
t' increases, conversely to what can be seen in Fig. 2. It 
is important to note that the horizontal axis in Fig. 4 is 
scaled by the Kondo temperature (for QDl) for t' = 0, 
estimated with the expression in the text above. In the 
inset to Fig. 4 we show that Tk for U' —Q and U' = U are 
equal. Indeed, as illustrated by the agreement between 
the entropy curves (for T < Tk) calculated with t' — 
[(black) O curve for U' = U, and (red) O curve for U' = 
0]. 
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FIG. 5: (color online) Magnetic moment ksTx as function of 
temperature for various values of t' . All parameters are the 
same as in Fig. 4. 

In relation to the magnetic susceptibility, a few dif- 
ferences are quite apparent when U' is turned on: (i) 
for the three smallest values of t' , in Fig. 5, the peak 
in the susceptibility has moved down to « Tk, and the 
maximum value of the magnetic moment is now lower 
(0.4 for finite [/', as compared to sa 0.45 for U' = 0). 
The decrease in the magnetic moment comes from the 
two extra states (two electrons in each dot) that arc now 
part of this manifold of states, and which have zero mag- 
netic moment; (ii) as t' increases, the peak moves up 
in temperature (up to T = 10^ Tk) and the maximum 
value of the magnetic moment further decreases, down 
to fJ-1mp/{9f^B)^ « 0.3 for t' = 0.01, (in) for the three 
smaller values of t', the splitting of the single plateau in 
the magnetic moment translates into a higher tempera- 



ture shoulder (around T/Tk = 10^), which is not present 
for [/' = 0. This shoulder becomes also the maximum 
value of the magnetic moment for t' = 0.01 (the largest 
t' results calculated). 

However, the most important differences coming from 
adding U' are related to the behavior of Tk and T^. 
First, the value of Tk being used for the scaling of the 
horizontal axis is that obtained when t' — 0. Note from 
Figs. 2 and 3, that Tk is very weakly dependent on t' 
when U' = 0. That is clearly not the case for U' = U, 
where it can be seen a strong variation of Tk and T^ with 
t'. In addition, more interestingly from an experimental 
point of view, it is clear that the ratio T^/Tk increases 
by a few orders of magnitude for U' = U. 
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FIG. 6: (color online) (a) Entropy Simp and (b) Magnetic 
moment kBTx as function of temperature for various ratios 
of U'/U (between and 1, see legend) for t' /U = 0.0026. 
To preserve the p-h symmetry for each curve we set Vg — 
— U/2 + U' . All the other parameters are the same as in 
Figs. 4 and 5. It is clear from the results that the Kondo 
temperature for the second stage of the TSK effect is greatly 
enhanced (by a few orders of magnitude) when U' increases. 

To analyze in more detail the increase of T'j^ with U' 
we show, in Fig. 6, results for entropy [panel (a)] and 
magnetic moment [panel (b)] which clearly indicate the 
strong increase in T^ (by a few orders of magnitude) 
when the ratio U'/U varies from to 1 for a fixed value 
of t' /U = 0.0026. This increase in T^ can be under- 
stood by estimating the Kondo temperature T^ for the 
second stage of the TSK, by using cq. 17, and explicitly 
calculating J' for finite U': 



J' = - 



U -U' - yJ{U -U')'^ + l&t''^ 



(18) 



indicating that, at least for small values of U'/U (where 
J' ~ 4i'^/(C/ — U'), there should be an exponential in- 
crease of T'^ with U'/U. This expectation is supported 
by the NRG results shown in Fig. 6. One can go one 
step further and, substitute eq. 18 into cq. 17, and use 
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FIG. 7: (color online) Fitting (dashed line) of the T'j^ values 
[red dots] extracted from Fig. 6, obtained by using eqs. 17 
and 18. The values obtained (from the fitting) for the two 
free parameters a and h are indicated in the figure. 



the equation thus obtained to fit the values of T^ that 
can be extracted from Fig. 6. This fitting is shown in 
Fig. 7. Note that the only free parameters in the fitting 
shown in Fig. 18 are 'a' and 'b' (values indicated in the 
figure). It can be clearly seen that fitting (at least up to 
U' /U ^ 0.7 is very good (note that the vertical axis is in 
a logarithmic scale). Indeed, this is one of the principal 
results of this paper. Given the recent advances in the 
use of floating interdot capacitors (see Ref. 32), values of 
U' /U for double dot systems have been steadily increas- 
ing, and may, in light of our results, offer the hope of 
observing the so far experimentally inaccessible energy 
scale T^. 



2. Zero-temperature case: QD occupation and conductance 

We first study how the occupation of the dots, (n^) = 
X]cr('^io-), is modified due to inter-dot Coulomb repul- 
sion. In Fig. 8, we show (ni) (top) and (712) (bottom) 
as function of Vg for various values of t'. For t' = 
[(black) O curve], as Vg sa one can observe a smooth 
increase of (rii), while (712) remains zero all the way down 
to Vg/U w —0.5, where (712) jumps to 1 and (ni) abruptly 
decreases. The discontinuous increase of (712) is due to 
the sharp peak in the local density of states (LDOS) p2 
of level €2, which is not broadened for t' = 0. On the 
other hand, (ni) increases smoothly, as ei is broadened 
by Fi. The value of Vg for which the first discontinuity 
of (77,2) occurs (expected to be at Vg « 0, when U' = 0) 
depends on the additional capacitive energy necessary 
to put electrons in different dots, due to the inter-dot 
Coulomb repulsion U'. 

As discussed by three of the current authors in a pre- 
vious work,^" this gate-voltage-dependent charge oscil- 
lation can be understood as a competition between the 



FIG. 8: (color online) Total occupation (ni) (top) and (712) 
(bottom) as function of gate voltage (Vg) for t' /U = (O), 
t'/U = 1.0 X 10"^ (D), t'/U = 2 X 10"^ (O) and t' /U = 1.0 x 
10~^ (A). Other parameters are t'/U = 5 x 10"* and V = 
3.2 X 10""^. The regions delimited by the vertical lines refer 
to QDl states, namely, empty (E), intermediate valence (IV), 
Kondo (K) and doubly occupied (DO) intermediate valence. 



Kondo and Intermediate Valence (IV) regimes: level £2, 
when t' = 0, acts as a dark state, whose charge occu- 
pation, which can only be an integer, adds a step func- 
tion of height U' to the gate potential Vg in QDl (de- 
pending on (712) being or 1). The value of Vg for 
which the transition of (772) between and 1 (and vice 
versa) occurs depends on Fi and reflects which many- 
body regime, Kondo or IV, better optimizes the energy 
of QDl. The dark state £2 acts like a switch between 
the two regimes and may have applications in quantum 



computation 



50,51 



For example, when the first disconti- 



nuity occurs, for Vg/U w —0.5, QDl is in a Kondo state, 
however, further decrease of Vg makes it more favorable 
for QDl to be in an IV state, which can be accomplished 
by charging QD2 (by exactly one electron), this, due to 
the capacitive coupling, increases the effective gate po- 
tential of QDl, discharging it, and bringing it back to the 
IV regime. As V^/f7 further decreases below —0.5, QDl 
starts to transition from an IV to a Kondo state; when 
Vg/U ~ — 1, again the IV state for QDl is more favorable; 
this regime can now be achieved by completely avoiding 
the Kondo regime through the discontinuous charging of 
QDl by one extra electron. Total charge is kept constant 
by discharging QD2 completely ((772) = 0). Finally, fur- 
ther decrease of V^, bringing it close to the p-h symmetric 
point [Vg = —{U/2 + U') = —3U/2], increases the charge 
in QDl to almost 2 electrons; at this point, the Kondo 
state in QDl is more favorable, and a new discharging 
and charging occurs of QDl and QD2, respectively, after 
which, each QD hosts one electron. These gate-voltage- 
dependent occupancy oscillations are clearly observed for 
all curves in Fig. 8, with the difference that, for finite t', 
there are no discontinuities. Indeed, for finite t', disconti- 



nuitics in the {ui) are no longer observed, we rather notice 
that the jumps are smoothened out as t' increases. This 
continuous charging of QD2 now results from the broad- 
ening r2 of the local bare (and also the many-body) level 
at QD2. For the non-interacting case it is easily shown 
that r2 ex t'^Ti. It is important to note the differences 
and similarities between this model and the one studied 
in Rcf. 50: in the latter, one has a two channel system, 
where the dark state, for small values of t^, acquires a 
finite broadening and smoothen out the discontinuities 
seen in the (rii) and (712) curves in Fig. 8. In the model 
being studied here, a similar process occurs: the dark 
state in QD2 acquires a broadening through its connec- 
tion t' to the single conduction channel through QDl. 
Nonetheless, through a comparison of Fig. 8 in this work 
with Fig. 4 in Rcf. 50 [panels (c) and (d)], one can see 
that qualitatively the results are very similar, indicat- 
ing that the basic processes determining the gate voltage 
dependent occupancy oscillations are the same. 

The strong gate-dependent variations in {rii) {i ~ 1, 2) 
(caused by the capacitive coupling) are expected to have 
dramatic influence in the conductance of the system, spe- 
cially in the Kondo regime. For [/' = 0, and in the special 
case where i' = 0, QDl will be in a Kondo state for all Vg 
values in the interval [—U, 0]. For finite U' , however, the 
charging and discharging of the QDs moves QDl from 
a Kondo to an IV regime, and vice versa. For the case 
of t' = [(black) O curve in Fig. 9], when Vg becomes 
negative, we observe a narrow Kondo plateau of height 
Go, and then the conductance drops suddenly to almost 
zero at Vg/U « —0.5, which is exactly the gate poten- 
tial value where QDl is discharged (sec Fig. 8). This 
drop signals the transitioning of QDl from the Kondo to 
the IV regime, as described above and in Ref. 50. The 
small peak (of height 0.3 Go) observed at Vg/U « — 1 
corresponds to the discontinuous jump of the renormal- 
ized level ei « ei + U'{n2), caused by the discontinuous 
change in (712) from 1 to 0, leading QDl from a low oc- 
cupancy ((ni) « 0.0) IV regime, to a high ((jii) w 2.0) 
IV regime, completely skipping the Kondo regime (with 
occupancy (ni) w 1.0). 

When the system reaches the vicinity of its p-h sym- 
metric point [Vg = -(C//2 + U') = -3[//2], there are 
exactly two electrons in the QDs (one in each) and the 
Kondo effect is fully reestablished in QDl (for t' = 0), 
thus the conductance reaches Go. Due to p-h symmetry, 
the other half of the curve can be explained in a similar 
way. The (black) O curve in Fig. 9 should be compared 
to the right-side panel of Fig. 3(b) in Ref. 50. 

For increasing t' there are mainly two differences in 
the charge and in the conductance: i) the abruptness 
of the charging/discharging and the discontinuities in 
the conductance are progressively smoothed out, and 
ii) the central peak in the conductance (around the p- 
h symmetric point) is fully suppressed. For example, for 
t' /U = 1 X 10~^, a rapid drop in the QD occupations is 
still clearly visible (as observed on the D (red) curve in 
Fig. 8), resulting in a rapid variation in the conductance. 
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FIG. 9: (color online) Zero-temperature conductance as a 
function of Vg for the same parameters as in Fig. 8. The 
occupancy discontinuities analyzed in Fig. 8 are associated 
to the discontinuities seen here for the conductance. See 
text for a full description of the variation of the peak struc- 
tures with increasing t': (a) t' /U = 0, (b) t' /U = 10~^, (c) 
t'/U = 5 X 10"^ and (d) t'/U = 10"^ 



which now moves down to a lower value of Vg when com- 
pared with the f' = curve (see D curve in Fig. 9). Notice 



that the peak at Vg 



-U is enhanced, as t' increases. 



as now the value of Vg in the IV regime corresponds to a 
higher conductance. For even larger t' there is a contin- 
uous enhancement of the conductance as Vg approaches 
—U, as clearly observed for t'/U = 1.0 x 10~^. The 
suppression of the conductance around the p-h symmet- 
ric point for small (but finite) t' results from destruc- 
tive Fano-like interference due to the Kondo resonance 
in QD2. For large t' , on the other hand, results from 
the formation of a local singlet due to an antiferromag- 
netic effective coupling (J = it'^/U) between the two 
electrons in the QDs, competing with the Kondo singlet 
formed between the electron in QDl and the conduction 
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FIG. 10: (color online) Conductance as function of V, for ^^G. 11: (color online) Conductance as function of tempera- 



C/2 = f/i, t'/U = 1.0 X IQ-^ for T/Tk = 7.31 x 10"!° (o), 
T/Tk = 1.08 X 10"* (D), T/Tk = 1.70 x 10"^ (O), T/T^ = 
1.66 X 10"^ (A), TITk = 6.49 x lO" (t>). 



ture for t'/C/ = 1.0 x 10"^ for VgjU = -1.5 [(blue) O] and 
Vg = —0.5 [(red) D], along the (blue and red) dashed vertical 
lines shown in Fig. 10. 



electrons. This is clearly associated to the TSK effect 
previously analyzed in this system. ^^ 

One should note that there is a striking difference in 
the effect of t' over the conductance for different re- 
gions of gate potential, and for different values of t'. For 
the smallest value studied (t'), as mentioned above, the 
Kondo peak at half-filling is immediately suppressed and 
this can be associated to the TSK effect. On the other 
hand, the peak around Vg = —U \% barely affected for a 
small t' . Indeed, as the TSK effect depends on the en- 
ergy scale J', one expects that it will be more effective at 
half-filling. Higher values of t' will then quickly modify 



the structures around Vg ~ —U/2: the peak at Vg 



-[/ 



is quickly enhanced and shifted to higher gate voltage 
values, while the discontinuity located (for t' = 0) at 
Vg = —U/2 moves to lower gate potential values, extend- 
ing the Kondo plateau, until the two structures merge 
and form a single Kondo peak centered at Vg = —U/2. 
This occurs because now level £2 acquires a finite width 
and therefore the discharging of QDl (and simultaneous 
charging of QD2) is not abrupt anymore and the system 
smoothly proceeds from the IV to the Kondo regime. 



3. Finite-temperature case 

Now we turn our attention to the effect of temperature 
in the transport properties of the system described in 
Fig. 1. In Fig. 10, we show the conductance for U' = U, 
t'/U = 1.0 X 10~^, and various values of temperature. 
The curves in this figure should be compared to the cor- 
responding zero temperature curves in Fig. 9 [(red) D 
curve]. Starting with T/Tk = 1-7 x 10"^ [(black) A] we 
notice a very similar shape when compared to the T ~ 
results [(red) D] in Fig. 9, except that the two symmet- 



ric Kondo plateaus and the IV peaks at Vg/U sa — 1 and 
Vg/U « —2 are slightly suppressed, which results from 
the small (but finite) temperature. Results for smaller 
temperatures (not shown) interpolate between T = and 
T/Tk = 1.7 X 10-3. For T/Tk = 1.66 x 10"! [(red) V] we 
observe a clear suppression of the Kondo plateau while, 
the IV peaks do not differ much from those for lower 
temperatures. This is readily understood since the char- 
acteristic energy of the Kondo state Tk is much smaller 
than the energy scale associated to the IV regime, of or- 
der Fi (for details, see Fig. 8 in Ref. 50). On the other 
hand, it is interesting to notice the emergence of a sharp 
small peak at the p-h symmetric point for a higher tem- 
perature such as T/Tk = 1.66 x 10" ^ [(red) V]. This 
peak results from a revival of the Kondo peak observed 
for t' ~ Q in the T = case, see the corresponding curve 
in Fig. 9 [(black) O]. This revival of the Kondo peak re- 
sults from the progressive restoration of the Kondo sin- 
glet state in QDl as the local singlet is suppressed when 
T^ <T < Tk- For T > Tk, as in the (green) > curve, 
we observe a suppression of the p-h Kondo peak, similar 
to a splitting. 

In order to show this effect in more detail, in Fig 11 
we plot the temperature variation of the conductance for 
the same parameters as in Fig. 10, for two different val- 
ues of Vg. For Vg/U = — 1.5 [(blue) O curve], the con- 
ductance vanishes when T —>■ 0, increases for 10"^ < 
T/Tk ^ 10"^, and is suppressed again for T > Tk- Note 
that the initial enhancement of the conductance as the 
temperature increases results from the destruction of the 
Kondo resonance in the QD2, suppressing the destruc- 
tive interference between the two paths. Conversely, as 
the temperature keeps increasing, and T exceeds Tk , the 
conductance is suppressed due to the destruction of the 
Kondo resonance in QDl. Finally, it will rise again and 



IV. CONCLUDING REMARKS 




FIG. 12: (color online) Conductance as function of Vg and 
t' /U = 0.01 for various temperatures in the U2 ~ Ui, t' /U = 
1.0 X 10-2 and T/Tk = 7.31 x 10"^° (O), T/Tk = 1.08 x 10""' 
(O), T/Tk = 1.70 X 10"=* (A), T/Tk = 1.66 x 10"^ (V), 
T/Tk = 2.59 x 10° (D), and T/Tk = 6.49 x lO" (t>). 



reach a maximum at approxim.ately T/Tk ^ 10^ ~ Fi 
due to charge fluctuations in QDl. For Vg = —0.5 [(red) 
D curve], i.e., at the edge of the right plateau of Fig. 10], 
the conductance starts from Go = 2e^//i as T — > 0, and is 
suppressed for 10"^ T^ < T < 10 T^ due to the destruc- 
tion of the Kondo effect in QDl. This is the behavior of 
the regular Kondo effect in a single QD. It is interesting 
to point out that once the temperature is above T^, the 
conductance at both Vg values behave alike. 



In Fig. 12, similarly to Fig. 10, we show G vs. K; 



for different temperatures, but now for i' = 1 x 10" , 
the largest t' value used in Fig. 9. These curves should 
be compared to the zero-tcmpcrature (blue) A curve in 
Fig. 9. For the case of T/Tk ^ 10-^°, the curve has 
the same shape as the corresponding curve in Fig. 9. We 
notice, however, that by gradually increasing the temper- 
ature, both Kondo plateaus become valleys surrounded 
by two Coulomb blockade peaks, while the gap at the 
p-h symmetric point remains almost unchanged. This is 
in consistent with what one would expect of the charge 
transport at high temperature for two molecular levels. 



We have studied a strongly interacting double dot sys- 
tem arranged in a T-shape configuration, where the dots 
are coupled via a tunneling matrix element and also 
strong inter-dot Coulomb repulsion. We have presented 
a detailed study of the effect of the inter-dot Coulomb 
interaction on the Kondo physics of the system. Our 
numerical analysis reveals interesting crossovers between 
different mixed valence and Kondo regimes that pro- 
duce dramatic changes in the conductance of the system. 
These crossovers can be tuned by varying the gate voltage 
of the QDs, and the charging/discharging processes pro- 
duce anomalous peaks in the conductance across the sys- 
tem, allowing a clear identification of the various regimes. 
We show that the inter-dot Coulomb repulsion not only 
preserves the TSK regime, but it also dramatically in- 
creases the lowest energy scale of the TSK effect. Indeed, 
T^ increases exponentially with the inter-dot Coulomb 
interaction, as shown in panels 9a) and (b) of Fig. 6. We 
stress that this enhancement may allow the experimen- 
tal observation of this so far elusive effect. Finally, by 
fixing the gate voltage at the p-h symmetric point and 
raising the temperature, the conductance shows clearly 
the crossover from TSK to Kondo and then to mixed va- 
lence regime (see Fig. 11). By contrast, away from the 
p-h symmetry only the regular Kondo regime is observed. 
For temperatures larger than Tk (and large enough val- 
ues of the coupling t' between the QDs) the conductance 
of the system as a function of the gate voltage possesses a 
four-peak structure, showing clearly the Coulomb block- 
ade regime of the molecular orbitals of the system. Fi- 
nally, we believe that our results will motivate future 
experimental measurements. 
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